Josephson and spontaneous currents at the interface between two <i-wave 
superconductors with transport current in the banks 
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A stationary Josephson effect in the ballistic contact of two 
d-wave superconductors with different axes orientation and 
with tangential transport current in the banks is considered 
theoretically. We study the influence of the transport current 
on the current-phase dependence for the Josephson and tan- 
gential currents at the interface. It is demonstrated that the 
spontaneous surface current at the interface depends on the 
transport current in the banks due to the interference of the 
angle-dependent current-carrying condensate wave functions 
of the two superconductors. 



I. INTRODUCTION 

It was shown that in the ground state of the con- 
tact of two d-wave superconductors with different axes 
orientation there is the tangential to the boundary cur- 
rent [1]- [8]. In the particularly interesting case of 7r/4- 
misorientation the ground state is two-fold degenerate: 
there are the tangential currents in opposite directions 
at 4> — ±7r/2 in the absence of Josephson current. The 
probabilities to find the contact in one of the two states 
are equal and the corresponding tangential current is re- 
ferred to as the spontaneous one. It was proposed to use 
such two-state quantum systems for the quantum com- 
putation [9]- [11]. It is of interest to study the possibility 
to control this system by the external transport current, 
which is the motivation of the present work. 

In the described problem of the Josephson contact of 
two d-wave superconductors with transport current in the 
banks the resulting tangential current is not a sum of the 
spontaneous current and the transport one. In the paper 
[12] we have studied simpler case of the contact of two 
s-wave superconductors with the transport current flow- 
ing in the banks. It was shown that the presence of the 
magnetic field [13]- [16], of the transport superconduct- 
ing current [12], or of the current in normal layer [17]- 
[18] in a mesoscopic Josephson junction can significantly 
influence current-phase characteristics, current distribu- 
tion etc. 

In the present problem the Josephson current is de- 
fined by the interference of the angle-dependent conden- 
sate wave functions of the two superconductors. There 
are two factors of anisotropy which define the angle de- 
pendence of the order parameter: the pairing anisotropy 



and the transport current. Thus, it is natural to expect 
that the resulting interference current (which has both 
normal and tangential components) is parametrized by 
the external phase difference (j> and by the value of the 
transport current (or by the superfluid velocity v s ). The 
presence of these two controlling parameters can be use- 
ful in the applications of Josephson junctions of high-T c 
superconductors. 

In Sec. 2 we derive basic equations to describe the bal- 
listic planar Josephson junction of two differently orien- 
tated d-wave superconductors with homogeneous current 
in the banks. These equations are solved analytically in 
Sec. 3. Then we study in Sec. 4 the influence of the trans- 
port current on the Josephson current and vice versa at 
the interface. In Appendix the order parameter and the 
current density in the homogeneous situation are consid- 
ered. 



II. MODEL AND BASIC EQUATIONS 

We consider a model of the Josephson junction as an 
ideal plane between two singlet (particularly, d-wave) su- 
perconductors with different axes orientation (see Fig.l). 
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FIG. 1. Geometry of the contact of two superconductors 
with different axes orientation and different transport currents 
(superfluid velocities v 3 -l,r) in the banks. 

The pair breaking and the scattering at the junction as 
well as the electron scattering in the bulk of metals are 
ignored. We did not take into account the possibility of 
the generation of a subdominant order parameter, which 
results in decreasing of the amplitude of the current [7]. 
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The c-axis of both superconductors is parallel to the in- 
terface. The direction of c-axis is chosen as z-axis. The 
a and b axes are situated in rry-plane. In the banks of 
the contact homogeneous current flows with a supercon- 
ducting velocity v s . We consider superfluid velocity v s in 
left (L) and right (R) superconducting half-spaces being 
parallel to each other v s l||v s /j and to the boundary; we 
choose y-axis along v s and x-axis perpendicular to the 
boundary; x = is the boundary plane. 

We describe the coherent current state in the supercon- 
ducting ballistic structure in the quasiclassical approxi- 
mation by the Eilenberger equation [19], [20] 



vf-k-G 
or 



u?3 + A,G 



0. 



(1) 



u n + ip F v Sj u„ = 7rT(2n + 1) arc the Mat- 

9 I 

f + -g 
A 
A* 

equation (1) should be supplemented by the equation for 
the order parameter (the self-consistency equation): 



where u 

subara frequencies, G = G w (v F ,r) = 
energy integrated Green function, A = 



is the 
. The 



A(v F ,r) = 7dVoT]T (V(v F y F )f(v' F ,r)) v , F 



(2) 



No is the density of states at the Fermi level and (■■■} v 
is the averaging over directions of vj?; V(\r F ,v' F ) is a 
pairing attractive potential. For the bulk d-wave super- 
conductor it is usually assumed A(i?) = Ao(T, v s ) cos 2i?, 
V(v F ,v' F ) = Vd cos 2i9 cos 2$', where angle $ defines a 
direction of the velocity \r F . Solutions of Eqs. (l)-(2) 
must satisfy the conditions for the Green functions and 
gap function in the banks far from the interface: 



9 (Too) 



/(Too) 



A(T-oo) 



^L,R 

A(=foo) = A i)J jexp(±i0/2). 



(3) 

(4) 
(5) 

Here w Lifl = u n +ip F v s . L ^ R , il LfR = yJul >R + A| fl ; 4> is 

the phase difference between the left and right supercon- 
ductors, which parametrizes the Josephson current state. 
The angles xl,r define the orientation of the crystal axes 
a and b in left and right half-spaces (see Fig.l). The an- 
gle between the axes of the right and left superconductors 
(the misorientation angle) is 5\ = XR ~ XL- 

Provided we know the Green function G, we can cal- 
culate the current density: 



j(r) = -2TTieN Q Tj2 (v F g(vF, r)) VJ? . 



(6) 



For singlet superconductors it is usually assumed 
A(— Vf) = A(vf) and therefore we have: 



f + (u,-v F ) = f + {-u,v F )=r(u,v F ), (7) 

g(u,-v F ) = -g(-u,v F ) =g*(u,v F ). (8) 
Making use of Eq. (8), Eq. (6) can be rewritten: 

j(r) = -jo^r J2 &f Im 5 (r)) VF , jo = ^ |e| N(0)v F T c . 



u>0 



(9) 



III. ANALYTICAL SOLUTION OF 
EILENBERGER EQUATION 

In this paper we consider the problem non-self- 
consistently: we assume the superconducting velocity v s 
being homogeneous and the order parameter A being 
constant in the two half-spaces: 



v s (r) 
A(r) 



v s;L , x < 
v s;i? , x > ' 
A L exp(i(f>/2), x<0 
A fl exp(-i0/2), x > 



(10) 



As it was shown in the papers [7], the self-consistent 
consideration of Josephson junction of d-w&ve supercon- 
ductors does not qualitatively differ from the non-self- 
consistent one. In the paper [7] the authors compare nu- 
merically the self-consistent solution with the non-self- 
consistent one. The self-consistency of the solution al- 
lows to take into account the suppression of the order 
parameter at the interface; the major effect of this is the 
reduction of the current [7]. 

Eq. (1) together with Eqs. (3-5, 10) yields for the left 
and right superconductors: 



WL,fl . n ( 2\x\ 
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/L,fl(x) = ^e- Si s n ^/ 2 



(11) 



(12) 
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where r\ = sign(u x ). Making use of the continuity con- 
dition, we obtain the expression for (/-function at the in- 
terface: 

_ SIlUr + £I r uj l - irjA L A R sm<j> 
Ql^r + ^luj r + AlA r cos <p 

Eqs. (9) and (13) allow us to calculate Josephson cur- 
rent jj = j x (x = 0) and the tangential current j y (x = 0) 
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at the interface. We underline that these equations are 
valid to describe the current at the interface of two sin- 
glet superconductors with different axes orientation and 
with different transport currents in the banks. The con- 
tact of conventional superconductors was considered in 
[12] and in the present paper we study the contact of 
d-wave superconductors, for which the order parameter 
is A L . R {d) = A (T,v s;Lii? )cos2( l ? - X l,r). The pre- 
sented here consideration can be also used to consider 
the contact of g- wave superconductors or s-wave/d-wave 
contact, etc. 

As we restrict ourselves to the non-self-consistent 
model we should calculate the order parameter Ao = 
A (T, v s ) in the bulk ci-wave superconductor. That is 
the subject of the Appendix. 

In the particular case considered below in detail we 
have v s; z, = w S ;R = v s and denote uj = u n + ipFV S} 



A\ R ; in this case we obtain 



5(0) 



Lu(flL + Qr) - ^AlA^ sin< 



A l Ar COS (j) 



(14) 
0) in this 



In the absence of the transport current (v s 
expression: uj = uj n [7] . 

We should also clarify the sign of the square root in 
^l,r- To make the solution (11) convergent, we should 
require RcQl,r > 0, which fixes the sign of the square 
root in Q.l,r to be sign(w • Pfv S :L,r.)- Moreover, this 
requirement, as it can be shown, provides the supple- 
mentary condition on Reg: sign(Reg) = sign(w). 



IV. INFLUENCE OF THE TRANSPORT 
CURRENT ON THE JOSEPHSON AND 
SPONTANEOUS CURRENTS AT THE 
INTERFACE 

Further we study the Josephson contact for the definite 
case: v sL = v s .r = v s and xl — and \R = i"/4. 

For small values of v s (in the linear \apFV s /T c approxi- 
mation) we can state the following approximate relations 
(which are valid for the values of (f> in the vicinity of 
±tt/2): 

j./(-v s , 4>) ~ jj(v s , <f>), j y {-v s ,<P) - -is/(v s , -<t>), 
and for difference Sj = j(v s ) — j(v s = 0) : 

5jj{~4>) ^ Sjj(<f>), 5j v {~(t>) ~ 5j v (4>), 
while at v s = 

3j{-4>) = -3j{4>), jy(-4>) = -3y{4>)- 

In the linear approximation the shift current 5j y is an 
even function of <p m contrast to j y (v s = 0). For the 
spontaneous current (at (f> = ±7r/2) the shift currents 
5j y are equal: 



j v (<f) = ±tt/2) =j s + Sj y , 
js(-n/2) = -js{*/2), 6j y (-7T/2) 



(15) 



8j y (n/2). 



In nonlinear consideration these shift currents are differ- 
ent for the two cases and it is discussed below. 

At Figs. 2-3 we plot the normal (Josephson) and tan- 
gential components of the current densities at the inter- 
face plane as functions of the phase difference 4> at low 
temperature. 
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FIG. 2. Josephson current density through the interface 
jj versus phase (p ( X l = 0, \R = f/4, T = 0.1T C ); 
Aoo = A (T = 0,v s = 0) = 2.14T C . 
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FIG. 3. Tangential current density at the interface j y ver- 
sus phase 4> ( X l = 0, XR = t/4., T = 0.1T C ). 

In the absence of the transport current: (i) j is an odd 
function of (f>; (ii) the normal component of the current 
(Josephson current) is 7r-periodic; (iii) in the equilibrium 
state at <j> = ±tt/2: jj = 0, j y (±w/2) = j s = =F \js\- In 
the latter case the tangential current exists in the absence 
of the Josephson current; for that reason it is referred to 
as the spontaneous current. The presence of the trans- 
port current breaks (i)-(iii) symmetry relations. There is 
non-zero Josephson current at <j> = 0, it. How the trans- 
port current influences the spontaneous current (i.e. the 
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tangential current at = 7r /2 and = — tt/2) is shown 
at Fig. 4. 
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FIG. 4. Fangential current density at the interface j y for 
the two values of phase difference (spontaneous current) ver- 
sus superfluid velocity v s (\l = 0, \R = 7r /4, T = 0.1T C ). 

The shift of the two values of the current for small 
values of v s (in the linear in ppv s /T c approximation) is 
equal (see Eq. (15)); however at values v s ~ 0.2A o/pf 
the shift current (i.e. the difference j y {v s ) — js{v s = 0)) 
is of different signs for the two currents and of opposite 
to js directions. 

We also note the following relations for v s 7^ 0: 
(1) 3j{4> = I") = -jj(<t> = 0) 7^ (the presence of 
the transport current induces non-zero Josephson cur- 
rent in the absence of an external phase difference); (2) 
3 j (0 = ±f ) =0, $$(<j> = ±f ) > (the transport cur- 
rent does not change the values of equilibrium phase dif- 
ference, at = ±7t/2); (3) j y {<j> = tt) = j y (<f> = 0) ^ 0. 
The latter relation concerns the interesting phenomena, 
studied in [12]: for some values of phase difference (here 
in the vicinity of = 0, tt) the interference of the angle- 
dependent condensate wave functions results in the ap- 
pearance of the additional tangential current with the 
direction opposite to the transport current in the banks. 
We underline that the resulting tangential current is not 
a sum of the spontaneous current and the transport cur- 
rent [12]. Thus, the transport current drastically influ- 
ences both the tangential (spontaneous) and Josephson 
currents. 

We can write down explicitly simple expression for the 
current for temperatures close to the critical (so close 
that A Q ,p F v s <C T c ). From Eq. (14) we have: 

Im0(O)£2A L A fl [-i7-^sm<£+5^cos0+ (16) 

3 ( Pf v s ) 2 . A L A R . 

+77- — sin + 77 4 sin 20] . 

At xl =0 and xr — ? r/4 this results in the following: 



3024tt ' T4 



~ 3 . A 2 {vfVs? . , 

3 = 56^ jo if^r sm<l) - e y 



(17) 
(18) 

(19) 

(20) 



jj = -1.7-10-^o sin 20 



is = -6.6 • 10" 2 io 1-77T sin0. 



Here A = A (7> s ) and is defined by Eq. (25). In 
particular, at v s — this gives: 



(21) 
(22) 



We note that j = — ^ PF ^ js- It follows that the 
effect of transport current on "spontaneous" tangential 
current at T ~ T c is to reduce its value by a small shift. 
It is remarkable that the tangential to the boundary cur- 
rent contains only corrections of the second order on the 
parameter pfv s /T c 1 . If xl = and xr — $X 7r /4, the 
integration of the second term in Eq. (16) would give us 
the factor tt cos 2 $x — § , which is zero for 6\ = \ ; this 
term at 8\ = and = gives the homogeneous current 
density (Eq. (26)). 

The integration of the first term in Eq. (16) gives us 
the factor cos 2Sx for x-component of the current and 
sin2c>x for y-component. In the case of 5x = tt/4 this 
term gives only the tangential component. As the conse- 
quence js > jj (see Eqs. (21)-(22)). 

It was discussed above that the linear in pfv s /T c terms 
result in the uniform shift of js ■ We can see that nonlin- 
ear terms result in the shift of different sign, and in both 
cases of opposite to js direction (see Eq. (20)). This in 
part explains the nonmonotonic behavior of j y (see Fig. 
4). The fact, that the presence of the transport current 
significantly changes the tangential (spontaneous) cur- 
rents, can be proposed to be used for its control, which is 
important in view of their possible application for quan- 
tum computation [9]- [11]. 

V. CONCLUSION 

The influence of the transport current, which flows in 
the banks, on the stationary Josephson effect in the con- 
tact of two d-wave superconductors is studied. We have 



1 There is also a term with the factor p ^" s which is ne- 
glected here. This term results in the equal shifts of js for 
<t> = ±7r/2. 



4 



derived the equations, which allow general consideration 
of the contact of two singlet superconductors with dif- 
ferent axes orientation and with different transport cur- 
rents in the banks. Particularly, we have studied the pla- 
nar contact of two d-wave superconductors in the case of 
7r/4-misorientation with equal transport currents in the 
banks. It was demonstrated that the current-phase rela- 
tion drastically depends upon the value of the transport 
current. The ground state degeneracy in the absence of 
transport current (at <f> — ±7r/2) is removed at v s =/= 0. 
The dependence of the shift current (which is the differ- 
ence of the tangential current and the spontaneous one) 
on v s is shown to be non-linear. It is proposed to use the 
transport current for the control of qubits based on the 
contact of two d-wave superconductors. 

We acknowledge support from D-Wave Sys. Inc. (Van- 
couver). 
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VI. APPENDIX. ORDER PARAMETER IN THE 
HOMOGENEOUS CURRENT STATE 

In this section we study the homogeneous current state 
in the bulk d-wave superconductor (see also in [21]). We 
note, that the order parameter Ao is the function of tem- 
perature T, superfluid velocity v s , and the angle \ be- 
tween the crystallographic axis a and the direction of 
superfluid velocity v s . For that we should solve Eqs.(2) 
and (9) with g and / given by Eqs. (3) and (4): 



- = 2T 



ir/2 



d-dRc 



A 2 Q?)/Ag 

n 



jo 



T 

nT c 



tt/2 

*» o -i/2 



ddvF Im 



Here A = N V d , uj = u> n +ip F v s , ft = ^/lo 2 + A 2 , A(t?) = 
A (T,v s )cos2(tf- X ). 

For T = (substituting ttT^ w by the integral J du>) 
we obtain the equations for the order parameter A and 
the current density j: 



v s Pf 



AW 



/ v s p F \ 
(23) 



J_ 

jo 



1 



VsPF 



Air T n 
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(24) 
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JL|««|C0B*| 



V S PF 



In Eqs. (23) and (24) the integration is performed in 
the region where A(ij) 2 < (v s p F ) 2 for i? e (-f , f ). 

At Figs. 5-6 we plot the order parameter A (T, v s ) and 
the current density versus superfluid velocity v s for dif- 
ferent angles \ a t low temperature. For comparison we 
also plot the curves for the s-wave superconductor. 
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FIG. 5. Order parameter Ao(T, v s ) in the bulk d-wave su- 
perconductor vs. superfluid velocity v a for different angles \ 
between v s and a-axis (T = 0). 
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FIG. 6. Current density in the bulk d-wave superconductor 
vs. superfluid velocity v s for different angles \ (T = 0). 

The numerical analysis at low temperature shows that 



here A 00 = A (T = 0, v s = 0) = & c e- 2 / x , £ = Ae' 1 ' 2 , 



2 in spite of the strong anisotropy of the pairing potential, 
- tie! order parameter A , the critical velocity vf and the 
crj/ical current j c depend weakly on the angle % between 
v s and crystallographic a-axis (see Figs. 5-6 and in Rcf. 
[21]). Namely, the respective difference is maximal for 
X = and \ — 7r /4 and does not exceed 0.1. For small 
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values of supcrfiuid velocity, i.e. in the linear approxi- 
mation on the parameter v s pf/T c , both Ao and j do not 
depend on \- 

For a temperature close to T c = fiLo c e~ 2 l x = 0.47Aoo, 
where (5 = ^e c (C is the Euler constant), both the gap 
function A and current density j, do not depend upon 
angle X- 

^IgfK 1 -!)-^ 2 ' (25) 

3 _7((3)Alp F v s 

jo ~ 32^ T? T c ■ (2b) 

The temperature dependence of critical velocity vf fol- 
lows from Eq. (25): ^f- = y^J^l - ^. 
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